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THEORY OF THE JEE SYPHOH 
By B. Szczenlowskl 

SUMMARY 


A new approach to the theory of the mixing of two cvirrents in an 
injector is presented which deals with an incompressible ideal fluid. 

The equations of continuity and of motion of a single flow (assumed 
one dimensional), althoiagh well known, are demonstrated in order to 
facilitate the understanding of the modified theory. 

Without advancing any new hypotheses as to the physical nature of 
the mixing phencmena, the theory shows new potentialities in an appro- 
priate shaping of the form of the walls of the mixing zone so as to 
improve the jet-syphon efficiency beyond that heretofore theoretically 
predicted. This inprovement is shown in spite of the fact that no math- 
ematical optimum conditions theoretically exist. 

A few examples of ways to Improve jet-syphon efficiency are indi- 
cated. Application of the new theory to the case of a constant-pressure 
ejector thrust augmenter is presented in an appendix. 


INTRODUCTION 


The theory set forth in this report is an attempt to give a new 
viewpoint on the problem of the mixing of two currents in an injector 
(l.e., jet syphon), as ccm 5 >ared with the theories now existing and 
accepted. This new approach is intended to clear up some ambiguity in 
the presentation of the problem in certain textbooks and papers on 
hydraulics in idalch the inte 2 T>retation of the law of conservation of 
momentimi seems to be oversimplified, in that the external forces (more 
partlcvilarly, reaction of the channel wall) are disregarded in con- 
stnictlng the mathematical expression for the balance of forces and 
quantities of motion in the general case of flow with variable pressiire 
along the mixing zone. 

The actual phenomena occurring in the channel during the mixing are 
very complicated and cannot be exp3J.citly elucidated without the intro- 
duction of certain physical hypotheses concerning the behavior of the 
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microscopic particles of the liquid. It is not intended to give in this 
paper such a detailed physical explanation of these actml phenomena. A 
few general deductions as to the behavior of the current, its velocity, 
pressure, and so forth, may be drawn, however, from the general princi- 
ples of conservation (of mass, of ener^, and of momentum), if some addi- 
tional abstractive ass\m5>tions are made. This is done in the section 
"Introduction to Hew Theoiy of Jet Syphon." An outline of certain prin- 
ciples of hydr odynamics , namely, the equation of continuity and the equa- 
tion of motion, is presented first to facilitate the understanding of the 
modified theory. 

Ells work was done at the Ecole Polytechnique (iMlversity of Ifontreal) 
and has been made available to the National Advisory Committee for Aero- 
nautics for publication because of its general interest. 


SYMBOLS 


A 


a 

c 


^ot 


AE 


tot 


F 


g 

H 


n 


P 

r 

s,x,l 

ds 


cross-sectional area of channel, sq ft 
acceleration, ft/sec^ 
velocity, ft /sec 

total (kinetic plus pressure) energy of fluid, ft -lb 
total energy loss due to shock 

external mass force, lb 
terrestrial acceleration, ft/sec^ 
punp total head, ft 

Wt 

weight output ratio; n]_ = — — for primary stream and 
n2 = 1 - n^^ for secondary stream 
static pressure, Ib/sq ft 
radius, ft 
channel length, ft 

infinitesimal element of channel length 
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Th ■fchmat, lb 

(Th)^ "bhrust of single (primary) stream, lb 

•k 

V weight oiitput, Ib/secj w^ for primary and W£ for secondary 

stream 

X potential of mass forces, ft^/sec^ 

z elevation in gravitational field, ft 

7 specific weight, Ib/cu ft 

external propulsive efficiency 

(%x) external propulsive efficiency of a single stream 

injector efficiency 

p density, 7/g 

T time, sec 

V ^ 

(D angular velocity, sec“^ 

V 

Dimensionless coefficients: 



« = Pl/Po 
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. ^ appendix only) 

\[^ 

Siibscrlpts; 

0 atmospheric conditions 

1 initial cross section of mixing zone 

2 final cross section of mixing zone 

3 diffxiser exit 
Superscripts ; 

' primary stream 

” secondary stream 

EQUATION OF CONTINUITY OF A SINGLE FLOW 


The following theory is based on the assumption that the flow is 
one dimensional, that is, that there is only one geometrical dimension, 
the length of the flow path. This unique geometrical coordinate is, how- 
ever, considered as curvilinear; that is, the "axis" of the flow may be 
shaped arbitrarily in the three-dimensional space, provided the flew vein 
is sufficiently fine to consider the pressure and the velocity uniform in 
the whole arbitrarily chosen cross section of the vein. Such an approach 
is, of course, not entirely exact, but it has given fully satisfactory 
results in dealing with the thermodynamic problems of stationary nozzles 
and diffusers. 

The other assumption made herein is that the fluid is ideal, that 
is, nonvlscous and involving no friction at the channel walls. 

Consider a certain flow vein whose axis is arbitrarily shaped in 
space (fig. 1 ) . Let its cross section A be variable with the length s 
of the flow path and with time. Cutting this vein at any point whatsoever, 
perpendicularly to the axis, with two planes A and B which are at an infi- 
nitely small distance ds from each other and stq)poslng that in plane A 
the cross section is A; the pressure, P; the velocity, c; and the den- 
sity, p; one has in the B plane: 
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Cross section: 


Pressure: p 


A 




Velocity: 


Density: 


c + 


P + 


ds 


Ss 


bs 


ds 


ds 


ds 


where the cross section is measured In square feet; the pressure, in 
pounds per square foot; and the velocity, in feet per second; and the 
density is 


P = 7/e 

where 7, in pounds per cubic foot, is the specific vreight; and 
g = 32.17^ ft/sec^ is the terrestrial acceleration. 

During the infinitely small time interval dr the amount pAc dr 
of the mass of fluid enters the vein element ds throTi^ cross sec- 
tion A, while another amount of mass 

(p + |Ed.)(A + |aB)(<= + |d.]ar 

leaves the vein element through the other cross section B. 

On the other hand, during the same time interval dr, the mass of 
the vein sector ds has increased frcm pA ds to 

as 

One has, therefore, arrived, xmder the condition of the continuity of 
flow, at the equality 
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ds - pA dB =i 


pAc dT - ^ ds^^A + ^ ds^^c + ^ dsj dr 


as the given amount of mass cannot change. 

Neglecting the infinitely small quantities of the second and third 
order and introducing the definition of velocity 


there is obtained 



A(ap) + A(Acp) = 0 (1) 

OT OS 


which is the mathematical form of the principle of continuity of flow 
and is called the equation of continuity. 

Particular cases may be easily deduced from equation (l) . For 
instance, if the fluid is incan5>resslble, p = Constant, and therefore 


+ A(ac) = 0 
or os 


Another particular case occurs when the channel cross section is 
constant, that is, when A = Constant, and therefore 






= 0 


If, at the same time, the fluid is incompressible, one obtains 
Bc/Bs = Oj that is, the velocity c is only a function of time (i.e., 
at a given moment it is constant throu^ the whole channel) . 

A further possible particular case is tiie steady flow, which is 
Independent of time; that is. 


~(Acp) = ^(Acp) = 0 
OS ds 
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Acp = = Constant = 

g g 


( 2 ) 


If, in addition, the fluid is incoanpressible, then y - Constant = 7 q, 
arid one gets 


Ac = Constant = A^c^^ 


( 2 a) 


EQUATION OF MOTION FOR A SIHGIE FLOW 


The principle of conse 27 vation of momentum, (i.e., of quantity of 
motion) resulting frcm Newton's law which states that the force equals 
the product of mass and acceleration enables one to obtain the equation 
of motion. 

The forces acting on sector ds of the vein are as follows: 

(1) On the one hand, the res\iltant of pressure forces exerted on 
cross section A and on the lateral surfacej this resultant force simply 
equala the product of the pressure p and of the projection of cross 
section A and the vein lateral siarface on the plane of cross section B 
that is. 



(2) On the other hand, one has, similarly, 
and the resultant force in the positive direction of s is 

^ i ll 

Apart from that, the external force F, called mass force, may also 
act on the sector ds. This force is bound with the mass of fluid and 
it always possesses a certain potential X such that the conponent of 
mass force in any direction n, reckoned per the unit of mass, equals the 
partial derivative of X in this direction, that is, bK/^n. One, there 
fore, has in the direction of motion s; 


clFg = |X pA ds 
06 
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and the resultant force will be 


dFg - A^ds = ^pAds-A^ds 
os os os 

The acceleration, of the element ds in the direction of motion is 

a = ^ ^ ^.+ 8 c _ ^ Sc I 8 c 
dr Ss dr St 5s St 

and the mass of the element is pA dsj therefore, the equation of motion 
(also called Bernoulli's equation) will be 

pA ds ^ - A ^ ds = pA ds fc + 1^) 

Ss Ss \ Ss St/ 

or 


SX 1 ^ ^ ^ bc_ 

Ss P Ss Ss St 

For steady (l.e., stationary) flow there is obtained 


( 5 ) 


and, after integration. 




c dc 


X 



+ Constant 


(^) 


Among the mass forces, the following three may, above all others, 
actually occur: 

(a) Gravitational force 

(b) Centrifugal force 

(c) Shock force 


All these are "external” forces, in contrast with possible "internal" 
mass forces (e.g., viscosity). 
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The gravitational force is caused by the gravitational field of the 
earth and it actually oceiors in heavy, inconrpressible liquids (e.g., 
water) . Its presence is independent of the character of motion. 

The centrifugal force may occur only in rotational motion or in a 
combination of rotational and progressive motions. In purely progressive 
motion centrifugal force does not exist, as here the radial acceleration 
(of the Coriolis type) is nil. 

The appearance of shock forces in any actiial hydraulic or aerody- 
namic device is, in general, iinavoidable. When the fluid (either gas or 
liquid) stream strikes a rigid wall, either fixed or in motion, the shock 
force appears in the form of a reaction of the wall, and this is, of 
cotjrse, an external force, relative to the stream. The same reaction 
may occior if the liquid stream strikes another liquid stream, which will 
produce a reaction force external to the former liquid stream. 

As is seen, the general equation of motion (eq. (3) ) is independent 
of both cross section A and its changes in space and time, provided, 
as has already been pointed out, the flow vein is sufficiently thin. 

Equation ( 4 ) is integrable if the law of variation of density with 
pressvire is known and given. Insofan as gases are concerned, these varia- 
tions depend on pressure and temperat\ire, and therefore the character of 
thermodynamic evolution (i.e., process) must be determined. In the case 
of an inccHDpresslble fluid, 7 = Constant = 7^, equation (ij-) is immedi- 
ately integrable, giving 


X - 



= Constant 

2 


( 5 ) 


where the value of the constant is to be defined frcm the initial, final, 
or boundary conditions. 

The following examples of mass forces may be cited; 

(a) Gravitational; it acts in the direction of the z-axls (fig. 2 ) 
and is simply 


dF = -(pA ds)g 

where the negative sign means that the force is directed down. It follows 
that X = -gz; thus, for an incompressible fluid. 


z + 



Constant = z^ 
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(b) Centrifugal; it acts in the direction of the radius r and is 
proportional to the centrifugal acceleration 


a = 

where cd = 2 jtn /60 Is the angular velocity, n designating the number 
of revolutions per minute . Thus the mass force becomes 


dF = (pA ds)a 


that is. 


and therefore 


2 “ 1,8(X) ^ 


1,800 


= Constant 

2g 


This equation is known and applied to rotating, dynamically acting 
machines (turbines, punqps, etc.). 

(c) As for the shock force, the case of the exchange of shocks 
between the two flowing liquids, as occurs in a jet syphon, is the main 
object of study of this paper, and it is disciussed in detail in the fol- 
lowing section. 


IHTRODUCTIQR TO HEW THEORY OF JET SYPHON 

The mixing of the two cxirrents in an Injector (i.e., jet syphon) is 
a typical example of hydrodynamic shock provided the velocities of the 
two currents are initially different, that is, provided they are different 
at the first momfint of the meeting of the two currents. The following 
computation is restricted to the case of an Incompressible fluid in both 
currents. 

Let the exit of the internal channel I be placed at x = 0 (fig. 3)^ 
where its cross-sectional area is A^' and the velocity is Cj^' . At the 

same cross section the external channel E has a cross-sectional area A^” 
and the velocity of the flowing fluid is c^" ^ c-j_* . Both fluids are con- 
sidered to be incompressible and of the same specific weight 7 ^. 
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As was stated la the "Introduction," a few abstractive, hut logically 
admissible, assimaptlons must be made in order to make it possible to draw 
conclusions from the three general principles of conservation. These 
assumptions are as follows: 

(1) The two currents are considered to atomize, one into the other, 
and to intermingle so perfectly that the exchange of shocks is not only 
between the surfaces but also between the two whole masses. 

(2) Obviously, the "exfcenaal" shock force F" acting on the second- 
ary fluid, being merely the action of the primary fluid, is always in 
equilibrium with the "external" shock force F‘ acting on the primary 
fluid, which is merely the action of the secondary fluid; thus , F" = -F ‘ . 

(5) In any given cross section, for example. A, at distance x, the 
pressure throu^out the cross section, that is, the pressure p' of the 
primary fluid, equals the pressure p" of the secondary fluid, as it would 
be senseless to assume two different pressures in two fluids which are 
perfectly mingled one with the other. However, the corresponding veloci- 
ties c' and c" are different, until an possible shocks have been 
exchanged. This occiirs for the first time at cross section Ag at dis- 
tance I, Tdiere c* and c" become equal (e.g., 02). 

Of course, A = A* + A" and, according to the law of continuity, 
the wel^t outp^Ifcs will be w^^ and W2, defined as follows: 


^1 = ?'cAl’°l’ = 7o^’°' 

W2 = 7oAi"ci" = 7oA" c" 

According to the section "Equation of Motion for a Single Flow," the 
mass force is 


dF = 


Zo aM 

8 bx 


dx 


that is. 


— r ^ ^ = Constant 

7 oJ A 7 o 2 g 


provided that, in general, the cross-sectional area A is variable with 
the channel length x. 
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gets 


Applying this last equation consecutively to the two currents, one 


i. THl _ £l _ (c*) = 

,J A' ' 


.^2 


7o 2g 


Constant 


i . El . = constant 

7oJ A" ro 2g 


Differentiating and introducing A’ = and A" 

as the condition p" = p* = p, give 


_ ^^2 






— , as well 


F. =!i^^.!idc' 


y c’ 
'o 


g 


F" » is * ^ 
ro =" « 


As, however, P” = -F', this reduces to 

A.fe+:|)dp + ^dc' +^dc» 
7o\c^ c'7 ^ g g 


= 0 


( 6 ) 


One of the siinplest. solutions of this equation is obtained in the 
particular case of constant pressijre throughout the channel, that is, 
p = Constant = Pq» this case 


V 


V 


Wj^ dc’ + W£ dc" = 0 


Wj^c' + W 2 C" = Constant = + ^2'^l" 


The final velocity is here 



Wj^Cl' + W2C1" 
°2 = + wg 


( 7 ) 


( 8 ) 


V 
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Introducing the conditions of continuity, that is. 




A' = 


7oC‘ 


A" = 


’'0° 


A = A‘ + A" 


the following eacpression is obtained for the channel cross-sectional 
area: 


^ ^ (wici« + V2°l")^l - (^1^ - 

V’ [("^1°1* ^2^1) - '"I'"*] 

wherefrom, it may be deduced that A diminishes with the diminishing 
velocity c* and attains a minimum value when c' takes its final value 
given by equation ( 8 ) . The form of the channel may be deduced from equa- 
tion (9) only if an additional assunqjtlon is made as to the vaxiatlop of 
c* with the channel length. Of course, the simplest assumption will be 
that c* diminishes linearly with the increasing distance x from the 
duct entry, namely. 


c 


I 


c 


II 





-“l") 

(w^^ + 

Wg)! 

w^(ci' 

1 

0 


W 2 )Z 


where I is the total length of the duct. 


During the mixing of two currents of Initially different speeds, a 
certain loss of energy occurs- because of shocks which do not contribute 
to the exchange of momentum in the direction of flow. These shocks do 
certain work which finally degenerates into heat (increase of temperatijre ) , 
the molecTolar motion induced by these shocks becoming disordered. There- 
fore, the sxjrn of the kinetic and potential energies of the two fluids 
must grad\mlly diminish. This loss may be expressed as 
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AE = 




Thus, the total energy loss in the injector after the velocity becomes 
imiform eind equal in the two cvirrents is 

- a[(o')^ - ^/] - - (=x’f] - (»i * 


( 10 ) 


is 


The total energy of the fl\iid (kinetic energy plus pressure energy) 

®tot = + 3(°i" f * ti * "a) ^%o 


Thus, the efficiency of the injector will be 


AE. 


^ E. 


tot 


(12) 


tot 


provided the flow and the fluids are perfect, that is, provided no fric- 
tion or other losses are Involved. 

In the particiilar case of constant pressure throii^out the channel, 
as quoted above,’ p = Constant = p^, and 


..\2 


AE 


tot 


2g(w^ + V 2 ) 


E 


tot 


+ W2(=i"f 

Sg 


= 


(^1^1* + •*^2°l“)" 


^ (^1 + ^2) ^l(=l')^ + ^2(°l”)' 
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is. 


The other peirtlcular case is a tube of constant cross section; that 


A* + A” = ^ = A = Constant = A_, = 

•V ■ V O 


7 C 7 G 
'o o 




7 Ct • 7„Ct " 
1 'o 1 


In this case equation ( 6 ) becomes easily integrable, as one of the two 
unknown velocities, for example, c”, may be eliminated. One gets: 


c" = 


^ 2 °l'*^l”®* 


(w^Ci” + WgC'jc* - 




(^1 + ^2)^1* ^l" 
(^l°l” '^2pl* ) 



2 

Ti = -n •+ yo^l^2°l'°l'(°l' - °l”) 

2 ^1 / , \2 
S(wici + WgCj^' ) 

It may be easily demonstrated that the pressure Increases steadily 
with decreasing c', that is, with the channel length, attaining a maxi- 
mum value at the final moment, when c’ = c" = C2^ that is, at the end 
of the duct. At the same mcsnent, assuming p^^ = Pq, 
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AE. 


tot 


VlVg|yi(ci"f + - °l”)" 


2g(v^c^” + WgCi* 



It may be easily proved that the efficiency is in this case generally 
lower than that in the case of the constant-pressure injector as studied 
above. 


aDUDY CF GENERAL CASE 


Consider again the general case, that is, the case in which the duct 
cross section varies arbitrarily with its len^h. In this case the pres- 
sure will also vary alon^ the duct and its final value, at the duct end, 
will be, in general, different from the initial value. Thus, to get the 
desirable final value of the pressxare, the usual nozzle (or the usual 
diffviser, according to vdiether Pg is greater or less than the desired 

value) should be applied. 

It should be first remembered that the equation of motion applying 
to the general case has already been given as (eq. (6)) 




+ 



S 


I 



II 


= 0 


On the other hand, one has the following conditions of contlniilty of flow: 


A‘ = 


Wi 






A" = 


^2 

’' 0 °" 


A = A* + A”j 


(13) 



r 
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where the total cross-sectional area of the duct A is arhitr2irily 
shaped along the duct length (fig. ^): 


A = <p(x) 


Flitting equation (13) into equation ( 6 ) and eliminating A", the 
following expression is obtained: 


dp - 


W 2 


gA(A - A' )' 


dA. + 


1 

1 

1 

1 

<M 

dA' 

S 7 o 

(A - A ')2 

(A') 2 j 

A 


= 0 


which is not integrable until an additional hypothesis is made as to 
changes of c’ (i.e., of A* ) with the duct length. However, the pro- 
cedure may be modified as follows: Let the decision as to possible 

changes of c* with the duct length be deferred until the mathematical 
solution of equation ( 6 ) is found (it may be next assumed that, for 
instance, c' diminishes linearly with the Increasing duct length x), 
while it is assumed first that the following relation 



is arbitrarily given Instead of A = <p(x) . 
found that 


Denoting 



it is 


A" = A - A* 



c 


It 


W2C' 



dc" 


'^2^°')^' ^1] 

(c'lif - w^^^ 


dc' 


Substituting into equation ( 6 ) gives 


r 


M + -L 

St 




( 14 ) 
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whlctL is integraMe, provided the function __>lr(c' ) is given. There is 
obtained 



■where a is an aioxiliary variable used only for Integration. 

At the end of the duct, where the total equalization of velocities 
is first achieved, one should have: “ 


c 


ft 


c* 




= C2 




V-, + Wo 


^(P2 - Pi) 



Wg 2 

ja^Tlf* (a) + 

^1] 

do, 

f 

xi|f(a) - wj 

|2 

+(a) 


(15a) 


It might be easily shown that both the particiilar cases as quoted 
in the section "introduction to New Theory of Jet Syphon" (p = Constant 
6 ind A = Constant) may be deduced from equation (l^) or (15) . 


In designing an injector, the aim is to have the greatest possible 
efficiency for any given initial conditions Wg, °l")* 

The form of the function i|r(c’) should, therefore, be so chosen as to 
have the greatest possible "value of 


c 2 - 2 

Wg^ja^'Ca) + wj 

da 


^ [at(a) - ^ 

+(a) 

^1 

(r, * 1 ""2 ..f 


w^ + W^ 

V 1 / w^^ + Wg V 1 / 



(16) 


which means that in the numerator 
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J = 




(IT) 


the function >|f should he so chosen as to have J attaining a maximum. 


In order to solve this problem correctly, the calculus of variations 
should be applied. The following discussion is based on the theory as 
given in reference 1 (ch. XXXIV, pp. ^k^- 6 oo). 


The general problem of the calculus of variations is as follows: 
A form of the function f (x) is looked for such that the integral 


J = r ^F(x,y,y') dx 

is a maximum, where y = f (x) and y' = x^ and x^^ being two con- 

stant limits of integration. Such a maximum does not always exist. There 
are several conditions for the existence of a mavlTinmi (or minimum). The 
first of these is that the following equation (called Euler's equation) 
must be fulfilled: 




y" + 


8(y’) 

In the present specific 


y« + 

^y Sy« dx dy» 
case the expressions 


0 


F = 


y = f (x) = +(a) 

^ ^ ^ W2^(x^' + wi) 

y / v2 

y(xy - w^J 


may be substituted into Euler's equation. The final resiilt 
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V- + w 

t(a) = 


Is obviously trivial because it requires that c^’ c^'*. 

The conclusion is that no extremum exists in the present specific 
case and therefore no form of the function i|r may be found giving such 
an extremum. 3 ^here is no doubt, however, that a certain form of t may 
give a higher value of efficiency than seme other form. Thus, the 

alternative is to try several possible mathematical forms of i|f. These 
forms may be chosen so as to meet the simplest mathematical operations 
and yet to have a certain degree of generality by introducing a few 
arbitrary constants, Triilch could next be so chosen as to attain the hip- 
est value of Tj^ possible in the specific case. 

One of the possible forms may be, for Instance 


Wt Ct " + w^c, * . , ti 

^(a) = 1 A - (a)^^ 


(“i') 


°1 


where p ^ 0 . Another is 

■Ha) . Yoh * - “) 


where A2 = and p $ 1 . Furthermore, it may be possible that 

ilf(a) = a - ba 

where a and b are constants. 

One more possible form, which is as follows. 


i|f(a) = a + J (18) 

will be studied more thoroughly below, as it seems to be relatively very 
simple and to give satisfactory resiilts. It must be strongly en^hasized 
here that there is no proof whatsoever that the form chosen is the best 
possible. On the contrary, it is very probable that other forms of 
may be found, which, although probably more involved mathematically, will 
give still hi^er values of efficiency for the injector. 
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STUDY OF SPECIFIC FORM OF DUC 2 AS ASSUMED BY EOJJAIIOH (I8) 
Assxnaing the specific duct form (eg. (I8)) 

if(a) = a + I 

means that 



Therefore, 



In addition, it is ass\imed that 

Ag = 


where p. is an arbitrary constant whose value should he so chosen as to 
have the highest value of the efficiency under the given circum- 

stances. Thus; 


^ (vi + V2)gl'°l“ 

2 n(vi" + Vi') 


M.W2(ci* - Cl") 

|(wiCi" + W2Ci') 


C * f* 
Cl Cj^ 

[H 

Kci" + 

3°1*] 

+ 

1 
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b . (1^ - ^)("l + ) 

+ Vl' ) - ( "1 + 

Suljstituting the above values into equation (17) and assuming, for 
greater simplicity in writing, that 


w. 


= n. 


^1 ^2 ^ 


W2 

= n^ = 1 - n, 

Wi + Wg 2 1 


tt 

= CT (0 ^ 0 < 1 ) 


one gets, after integration, 

L (1 - n^Xl - - (1 - ni’o log, (■>! + ^-^) + 


_ 22 _ 


(«l') - ”lf(l - ”f 


( ^ 

2«2 j 


■^^(1 - >^f(l - ("I* + 1 - “l) 


+ 1)0 +1 - »i] l=Se (1 - ♦ 2 - log, n 

°l(l - °i)(l - ®)["l(3"i - l)» + C 1 - - 2)] - 


+ 1 - =1) log, (nj. + + 2(1 - + 1 - pj^) log, 


1 

7 


% = 


^/(°i' y 


»! + (1 - 


( 19 ) 


(20) 


Now the value of |x giving the highest value of t)^ should he 
found. Equation ( 20 ) is of the form 


Bo B.7C. 

^i = ^1 + IT ^ 

^ u 


<7 
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where a. 2 j and are constants. Thus 


dill ®-2 


2a? 


^.5 


dti' 


T)i 2ag _ 


The condition 


8*2 


8a.^ 


•f and 


^ = 0 will give ti =- 

®-2 . 

(^i)e = ^ - (^27^)- Thus a maxlmign of will exist only if 


d^Tii 


dfi' 


^ <0, that is, if a^ < 0 and, consequently, a^ > 0, because 


lig < 0 would have no physical sense. 


The correct mathematical analysis of the coefficients a^ and a^ 

is rather difficult j numerical testing leads, however, to the highly 
probable conclusion that a^ may never be negative, that is, nay 

never have a mathematical maximum, which means that it tends mathematically 
to an infinite positive value at p. = 0. One, therefore, should look for 
the highest value of in the whole range of physically possible values 

of p, which is not a mathematical maxi mum. 


It shoiild be noticed first that p cannot drop to zero, but it has 
a limit 


V-n = 


. KjJjgK” _ 


^1°1 


+ V2C3_« 


nio 


+ 1 - 


“1 


( 21 ) 


at which the duct cross-sectional area will becone infinitely large, 
except for initial and final values, if the trivial condition c^" = c^' 

is not fulfilled in addition. Thus, it must be assumed that 


and, therefore, the highest value of corresponds to this limit value 

of p. 
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Siibstituting |j.j as defined "by equation (21) Into equation (20), 
the limit value of the efficiency Is found: 


(^ 1)7 = 


2(1 - ^if loge (nj ; j +\ -- "n -) + °l(3^1 ~ + ^l(^ ~ ^l)^ 




It may he easily proved, for the sake of checking, that, for 0 ^ n^^ ^ 1, 

Cj^' ^ 0, and c^" ^ 0, which conditions are physically obvious, the value 

is always ) ^1. Substituting, namely, x = i , gives the 

' n^tr + 1 - n^ 

condition 


(x - 1 ) (3 - 2n3^)x - 1] 
2(1 - n]^)x 2 


^ logg X 


To prove that this condition is fulfilled for 0 < x ^ i , that is, 

1 ” n^ 

for 0 ^ a ^ 00 , it is sufficient to take the first two derivatives of the 
two sides of the inequality cited above. This inequality becomes an 
equality only for x = 1 . 

Of course, the extreme case as defined by equations (21) and (22) 
is not actually applicable^ some deviation of numerical values must, 
therefore, be applied and will become less than . In order, 

z 

however, to have con 5 )£trative figures on hand, the theoretical case of 
equation ( 22 ) will be studied numerically below, it being CQn 5 >ared with 
the value (^li^ s-s applied in the usual routine now (ref. 2 ), 

It may be easily proved that always 




It is sufficient, for instance, to substitute 
leads to the condition 


n^^cr + 1 - nj_ 


= X, which 
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(x - 1)(5X - 1) > ^ 


2x2 


TaMng the two first derivatives of this inequality, one arrives at the 


concliision that x ^ 1, which means that c-|_* 




This ass-umption was 


made at the very beginning and remains, of course, valid; therefore, the 


assumption that (ti^) = (%) 

If 

is, X = 1 when ( 


la proved, except for c-, 


= c 


1 f 


that 


KUMERICAL CCMPUDATION OF EQUATIONS (22) AND (25) 


The procedure for the numerical calculation of equations (22) and 
( 23 ) is as follows; Take first the routine formula for efficiency ( 


The numerical values of 


n^^, for different constant values of 


^ are given in table I and plotted against 

o. 


As is seen, ^ 


, in figure 5* 
passes throu^ a minimum value, the formula for 


which is easily obtainable by differentiation of equation ( 23 ) : 



(1 + af 


for 



a 

1 + a 


Numerical values of (ili) are given in table H. 

^mln 

Turning now to equation (22), it should be stated that for very 
small values of n^ the usual calculation based on form (22) is not 

sufficiently exact. It is, therefore, better to develop the logarithm 
as follows ; 


log^ 


1 


1 

(1 - a)n^ 


(1 - g)n2^ + i(l - a)^n^ + i(l - a)^n^^ + . . . 


which gives 
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( 11 ), 


o® + |<1 - o)(l + j + 0^)j^ + ^(1 - o^XJ- ^ ■»■ - g)^(l 4- » ■<■ 6g^)n^ ■!■ - b)**^! + 4j + 100 ^) 0 ^*'’ h 

+ (1 - o®)n^ 


(2i) 


The numerical values of (^i) given in table III and plotted 

against n^^, for_different constant values of a, in figure 6. 


DUCT FORM 


In order to deteiinine the duct form, an additional hypothesis as to 
changes of c' with the duct length must be made. As has already been 
pointed out in the section "Study of General Case," a possible assumption 
is that c' . diminishes linearly with the increasing duct length x (see 
fig. 4 ). Such an assumption seems to be the simplest and, at the same 
time, logical enough. Mathematically, it is equivalent to 


c' 



x+ 


( 25 ) 


In the particular case as defined by equation (l8) this leads to: 
ac* + b 

a(cg - + (aci‘ + 

’^o[(°2 - + °l'^] 

KVl" ^2°l*)[(^lV ^2^1*)^ - ’^2(''l* - 


(26) 

where p must, of course, be chosen a little larger than the limit value 
defined by equation (21), but close to it. In order to maintain the high- 
est possible efficiency. The length I will be evaluated according to 
practical considerations, as no theory may be advanced here in this con- 
nection. The aim is, of course, the thorough mixing with minimum losses 


’'o°l'°l" ■ ^"K°l'' + ^2°l’)^ - ^K°l" 


WgCi') - (wi + V2)cijxj 
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and the obtaining of ccsaplete unlformlzatlon of velocity at (or before) 
the final cross section. 


As for the cross-sectional area A, It may be assvmed In any form 
whatsoever. If a circular cross section Is assumed, the diameter will 



In the limit case p. = the results become, of course, trivial; 
that Is, 


A = A^ = A2 = Constant 


(wi + Wg) 


"1’ 


Cg = Cl" = Cl' 

The poirtlon I of the duct lying between cross sections A^^ and A2 

cannot be regarded. In general, as a whole unit, either in the case of an 
injector or in other cases. The reason Is that the "final" pressure Pg 

is. In general, not equa,l to p^. Heither does it equal the press\jre Pq 

of the surro-undlngs . In the case of an injector used for pun5>ing purposes, 
it may be assimied that p^ = Po^ 'tlie final pressure should be raised 

as high as possible, for exan^ile, to p^, by applying the feasible mini- 
mvmi c^ of the final velocity. This may be done only by adding a diffuser 
extending from A2 to A^. 

The other possible case Is a thrust producer. Here not only p^^ = p^ 
but also p^ = Pq should be assumed. 

In the general case, according to equation (16), 



and, according to equations (l5a) and (lT)j 

P 2 - Pi = P2 - Po = - =2^) 

wherefrom 

P2 - Pi = P2 - Po = § |ni[^i + (1 - ^l)'^](°i' f 
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In the additional diffuser (see fig. 7 ) the pressure is brought back to 
Pq, that is, p^ = Pg. Thus 


P5 - P2 (P2 " Po) ® 


— 2 . 

2 g 




But, according to Bernoulli's law, if friction and other losses are 
disregarded. 



Thus: 




(28) 


In the particTilar case as defined by equation (18) , the value 


CCt * 

^ p[l - (1 - cf)n^] 

should be entered into equation (27), and the value of as defined 

by equations (19) and (20) should be entered into equation (28) . 

The area of the final cross section of the diffuser is, obviously. 




7'o^5 




(29) 


The particular case of an ejector thrust augmenter with constant- 
pressure mixing zone is studied more thoroughly in the appendix. 


Ecole Polytechnique (tMiversity of Montreal), 

Montreal, Quebec, Canada, January 15 , 19 ^ 9 * 
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APPENDIX 

THEORY OF "CONSTANT-PRESSURE" EJECTOR THRUSC AU04ENTER 


In the computation developed in the main text the case of an injector 
has been treated almost exclusively, with the aim of obtaining the dynamic 
compression. The ultimate goal has been to get the largest possible total 
pressure (i.e., a sum of static pressure and velocity head). The case of 
an ejector thrust auigmenter is different in that the ultimate goal is to 
get the largest possible thrust, depending directly on momentum, under 
given initial and final conditions. This case is briefly investigated 
theoretically below. 


Althou^ the definition of efficiency is a little different in this 
case, it remains intimately related to the final momentum of the mixture 
leaving the ejector; therefore, all the general remarks concerning the 
ejector (mixing zone) throughout the paper remain valid. However, as the 
investigation of the most general case would be rather involved, it seems 
more practical to investigate first the particular case of the ejector, 
that is, of the mixing zone, at constant pressure. Therefore formulas (8) 
and (9), applying to the mixing zone, will remain valid. 


Obviously, if the given initial conditions Pq,Cq* for the primary 
stream and Pq,Cq" for the secondary stream are applied directly at the 
ejector exit, that is, if Pq_ = p^' = Pq"? *^i' “ °o*' °l" “ °o”^ 

and if, in addition, Pq ~ Pi ~ 5*2 atmospheric pressure, th\as no 

diffuser or nozzle at station 2), there is no thrust aiigmentatlon whatso- 
ever, because the momentum remains unchanged. It therefore follows that 
the given initial conditions should be first subject to appropriate changes 
if any thrust augmentation is to be obtained with a constant-pressure ejec- 
tor. These changes may be either a prior acceleration (by nozzles) or 
deceleration (by diffusers) before the two currents enter the ejector, 
with the condition p^’ = pj^" = p^^. Of course, in such a case a diffuser 

(or nozzle) must be applied downstream of the ejector, in order to estab- 
lish the desired final pressure p^ = p^. 


For the given initial conditions p^*,Cq* and Pq",Cq" and the 
final pressure p^, there exists a certain value of p^ = P 2 yielding a 

maximum of thrust, which may exceed that obtained with the primary stream 
alone and without ejector. Of course, these optimal values of p-j^ 

and Th depend on the mass-flow ratio, that is, on nn = = — . 

' ' -L + W- 


In order to define the problem further, assume the system as shown 
in figure 8> in which the boundary losses in the ejector and the exit 
diffuser are disregarded. 
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”1 

’'o 

H 


The given values are: 

atmospheric pressure, Ib/sq ft 

velocity of travel, ft /sec 


mass-flow ratio. 


w^ + Wg 


specific gravity, Ib/cu ft 


act\ml head produced by pump, ft 

The value of c, ' is to be found for ( Th) : thus values of p_ 
1 ' max 1 

and C2^" will thereby be defined. 

The following formulas follow: 


H 


■ (°o") ] - Sg[(“o') ■ (°o") j 


=1" = f f (Po - Pi) 


W^Ct ' + W„Ct" , V 

=2 ° ^ t Va = Vi' + - “i)”! 


" - ^1) 


Th = i(v^ + Wa)(oj - c„) 


g\ 1 d.!\ 3 

The (external) propulsive efficiency is 


c (Th) 

T1 = 

'ex 
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If the ejector were not applied, the primary stream alone would produce 
the thrust 



Introducing the dimensionless factors 

Pl/Po = « 




the following expressions are obtained: 
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As is seen, ^ is a function of fottr different factors, n^, A, 

and logically, they are limited as follows: 


0 ^ nj^ ^ 1 
0 = CTq = “ 

0 < X ^ 00 , 

0 < It ^ 00 

where it is understood that 

X < CT^' ^ ^X^ + + 1 

0 < It ^ X^ + cJq^ + 1 
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and therefore 


- (“1')^ “o' 1 i ° 

In order to determine the influence of these foiir factors, the for- 
mula for ^ should he examined first with respect to each of them sepa- 
rately. The condition d^dnj_ = 0 leads to 


^ n ^=0 


if It < 1 


(“ 1 ), = ° 

It may he shown that for n^ = 0 

(° i ' - /(“!')' - ^ 

except for ^i.e., A = 0 where 

_Q = 1. This is a maximum because, for n^ = 1, ^ = 1. Therefore, 

the secondary mass flow should he as large as possible. 

The condition d^ ~ ^ <ioes not lead to any extremal values, 

simply yielding A = 0, which does not impose any condition upon . 
It may he shown that, in general, for A = 0, 


^A=0 " 


= 1 




It may also he shown that, for p^^ = Po# 




‘^1' ‘^o 
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as well as that, for = 0, 


— 


1 

1 


(“1’) 

1 

CM 

l/M 

- A^ < 1 if A > 0 


Furthermore, ^ may become very hi^ for small values of (the 

efficiency being low, however), whereas it attains certain maximum valiie 
(greater than l) at a definite value of A. The picture is as shown in 
figure 9* 


Finally, it may be shown, in general, that ^>1 for 0<it<l 
and 0 < n^ < 1. 


If the ejector is not applied, then ir = 1, n^ = 1, 
and (T)ex) = * Therefore, the ratio of the increase 


+ ‘^o 


ciency by an ejector is 


= 1 , 

in effi- 



If the values of ^ot appHed without the ejector, are going 

to be changed (e.g., A will Increase) while remains the same, then 
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In the applications of a centrifugal puisp to the dynamic propulsion, 
the requirement of a reasonable efficiency Irgjoses a rather large mass 
flow while the head should be kept small. This results in huge tubing 
diameters . In order to improve this situation, an ejector thrust aiig- 
menter may be applied in such a manner as to keep the efficiency unchanged, 
while the pump head 1 \ is increased with subsequent reduction in mass 
flow, if the power supply remains unchanged. The problem may thus be 
defined as follows: For the given e* =1, increase \ and change at so 

as to maintain a staff iclently high value of thrust, while the velocity of 
travel remains unchanged. 

The condition =1 leads to 


(1 -*) 


4.(1 - ni)(oi^' + <Tof 


+ 2=1(01^' + t^o) (2 - 


(%' + “of - 3=1^) - '‘(1 - - °o^]° 


.wherefrom the functional dependence of at 


numerical values of 


"1^ 


and 


on A may be found for assumed 
The thrust will remain constant 


because the power supply lamains the same; that is. 
Therefore, 








+ CTr 



W^A^ 



w, \ 


= 1 

Let Pq = 2,120 pounds per square foot, 7 ^ = 62.5 pounds per cubic 
foot, Cq = 15 feet per second, Hq = 15 feet, otg = 1 and n^ = 0 . 25 . 
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Then cTq = O.52O, 7 \q = 0 . 662 , ct^q* 0 . 755 ^ ci^' = 5 ^. 5 ^ feet per sec- 
ond, Wg = and 

6 k 2 

1 _ rt = 0 - 01872 ?^ + 0 . 059 ^^A + 0 . 0767 A - 0.04485 

0.458 - 0,25A^ 


O'!* + (1 - jt) + 0.1024 

“ 0.755 . 

The reduction of the exhaust tube will he In proportion 

6 =,/ 5 ^ 

where 


r 

Some calcxilated values of 1 - « and 6 for given values of 
axe presented In the following table: 



Wi/ 

Wn 

■^0 


/%' 

8 




0.5 

0.6 

0.7 

0.8 

0.0542 

0.090 

0.169 

0.269 

0.876 

0.750 

0.626 

0.5475 

1.086 

1.210 

1.540 

1.472 

0.898 

0.776 

0.684 

0.610 



For the same assimiptlons but with n^^ = 0 . 2 , 


h p 

T 0.0U24 a° + 0.0508X + 0.07794 a - 0.04485 

X - It = — ■ " 

0.458 - 0 . 2 A^ 
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for the following veilues of 




0.438 

0.5 



0.8 

mm 


E9 

B 

Esa 

1 - rt 


0.0244 

0.0712 

0.1289 

0.2002 

0.2896 

0.400 

0.539 

0.713 

0.939 

"iZ-a, 

■ 

0.876 

0.750 

0.626 

0.3475 

0 . 48 ^ 

0.438 

0.398 

0.3648 

0.3368 


m 

1.078 

1.196 

1.312 

1.429 

1-547 

1.668 

1.795 

1.930 

2.080 

s 

1 

0.902 

0.782 

0.690 

0.619 

0.559 

0.5125 

0.471 

0.435 

o. 4 o® 


Finally, for = O.if, 

^ 0 «Q 599 X^ + o.o689^^ + o.o6l\^ - 0.0W3 

O.lj -38 - O.lt-A^ 


for the following values of 



0.438 

0.5 

0.6 

|Q|R||||^ 

0.8 

0.85 

1 - rt 

0 

0.0436 

0.1489 

0.3302 

0.667 

0.955 


1 

0.876 

0.730 

0.626 

0.5475 

0.5150 


1 

1.093 

1.256 

1.447 

1.704 

1.877 

6 

1 

0.895 

0.7625 

0.658 

0.566 

0.524 


The above numerical results are presented graphically in fig- 
ures 10(a) and 10(b). 
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TABLE I 

VALUES OP (ti^) 


\ 

(ll)r 


0 

mm 

1 

0.2 

0.3 


m 

0.6 

0.7 

0.8 

0.9 

1.0 

0 

0 

0.1 

0.2 

0.5 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.000 

0.1 

1.000 

0.551 

0.5768 

O.Mf 6 

0.5215 

0.599 

0.678 

0.758 

0.858 

0.919 

1.000 

0.2 

1.000 

0.576 

0.5585 

0.590 

0.6375 

0.6925 

0.750 

0.811 

0.874 

0.956 

1.000 

0.3 

1,000 

0.756 

0.712 

0.716 

0.741 

0.775 

0.815 

0.858 

0.904 

0.952 

1.000 

n k 

V »~r 

1 rv\pi 

rJ-1 « 

W ■ WW 1 ^ 

n fto=: 


0,826 

n RkliBi 

w *w • 

Q.870 

n.Aoo 

— - 

— m 

1,000 

0.5 

1.000 

0.931 

0.900 

0.890 

0.891 

0.900 

0.915 

0.9325 

0.953 

0.975 

1.000 

0.6 

1.000 

0.966 

0.9^75 

0.959 

0,938 

0.9415 

0.949 

0.958 

0.970 

0.985 

1.000 

0.7 

1,000 

0.985 

0.975 

0.970 

0.969 

0.970 

0.975 

0.977 

0.984 

0.9915 

1.000 

0.8 

1.000 

0.995 

0.9915 

0.9885 

0.968 

0.988 

0.9885 

0.990 

0.993 

0.996 

1.000 

0.9 

1.000 

0.9985 

0.998 

0.9975 

0.9972 

0.997 

0.9975 

0.998 

0.9985 

0.999 

1.000 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1,000 

1.000 

1.000 

1.000 

1.000 

1.000 
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TABLE III 
VALUES OP ( TJi ) 

It 




\ 

0 

0.1 0.2 0.5 0.4 0.5 0.6 0.7 

0.8 

0.9 

1.00 

0 

0 . 666 ( 6 ) 

0,6840 0.7050 0.7255 0.7460 0.7720 0.8015 0.8555 

0.8770 

0.9270 

1.0000 

fil 

1 _(V>nn 

o.Tion n.'TQTV) 0-7^ D-T^Tp; n.Tfip^ n.fl46o 

o.flfl4n 

O.QPPO 

T.onnn 







0.2 

1.0000 

0.7950 0.7726 0.7768 0.7890 0.8092 0.8550 0.8650 

0.9008 

n o)i)iii 

V/ » ^ 1' 1 '1 

1.0000 

0.5 

1.0000 

0.8686 0.8575 0.8250 0.8515 0.8470 0.8655 0.8900 

0.9200 

0.9565 

1.0000 

0.4 

1.0000 

0.9170 0.8880 0.8760 0.8760 0.8840 0.8970 0.9164 

0.9400 

0.9685 

1.0000 

0.5 

1.0000 

0.9520 0.9285 0.9180 0.9170 0.9214 0.9290 0.9412 

0.9560 

0.9765 

1.0000 

0.6 

1.0000 

0.9740 0.9600 0.9520 0.9495 0.9500 0.9550 0.9622 

0.9720 

0.9850 

1.0000 

0.7 

1.0000 

0.9770 0.9790 0.9750 0.9725 0.9726 0.9750 0.9790 

0.9840 

0.9915 

1.0000 

0.8 

1.0000 

0.9950 0.9910 0.9894 0.9880 0.9880 0.9890 0.9910 

0.9955 

0.9965 

1.0000 

0.9 

1.0000 

0.9990 0.9975 0.9960 0,9955 0.9960 0.9965 0.9980 

0.9985 

0.9990 

1.0000 

1.0 

1.0000 

1.0000 1.0000 1.0000 1.0000 1.0000 1,0000 1.0000 

1.0000 

1.0000 

1.0000 
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Figure 5 *“ Graphical representation of injector efficiency as con^juted 

with conventional formula. 
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Figure 6 .- A case of injector efficiency as con^iuted according to new 

theory. 






Figure 7.- Illustrative diagram of e .iec tor thrust augmenter. 



Figure 8.- Illustrative diagram of constant-pressure ejector- thmst 

augmenter . 


NACA TN 5385 


47 



Figure 9 «- Graph for efficiency and thrust of constant-pressure ejector 

thrust augmenter. 
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(a) Variation of a-^ 6, and wltli jt. 

Figure 10.- Dimensionless cliaracterlstics of constant-pressure ejector 

thrust augmenter. 




(To) Variation of and 5 with s. 

Figure 10.- Concluded. 
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